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Abstract 

> 

' We find the spectrum P{uj)duj of the gravitational wave background pro- 

duced in the early universe in string theory. We work in the framework of 
the recently discussed String Driven Cosmology, whose scale factors are com- 
Q^ I puted with the low-energy effective string equations as well as selfconsistent 

^ ■ solutions of General Relativity with a gas of strings as source. The scale fac- 

tor evolution is described by an early string driven inflationary stage with an 
instantaneous transition to a radiation dominated stage and successive mat- 
ter dominated stage. This is a string cosmology expanding evolution always 
. running on positive proper cosmic time. 

The spectrum of the generated gravitons is computed in the framework 
of Quantum Field Theory and in the appropiated effective String Cosmology 
context. We study and show explicitly the effect of the dilaton field, charac- 
5^ I teristic to this kind of cosmologies. We compute the spectrum for the same 

evolution description with three differents approachs. Some features of grav- 
itational wave spectra, as peaks and asymptotic behaviours, are found direct 
consequences of the dilaton involved in the string low energy effective action 
and not only of the scale factor evolution. 

We make use of a careful treatment of the scale factor evolution and in- 
volved transitions. This allows us to compute a full prediction on the power 
spectrum of gravitational waves without any free-parameters. A comparative 
analysis of different treatments, solutions and compatibility with observational 
bounds or detection perspectives is made. 
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1 Introduction and Results 



In the latest years, an impressive scientific effort has been made on the subject 
of production of rehc gravitational waves. The spectrum of relic gravitational waves 
constitutes an important source of information about the conditions and dynamics of 
the early Universe. From the pioneering works relating the very early Universe with 
the spectrum of metric fluctuations ([|1|, 0,0) a lot of work has followed computing 
the gravitational waves background in different evolution models (see for instance 
[i]5 0)05 0]'llll'[i])- "^^^ ^^^^'^ gravitational wave background is generated as a con- 
sequence of the amplification by the Universe expansion of the graviton fields, which 
decouples from matter at very early stages. Therefore, their spectrum is defined by 
conditions in the very early epoch and remains practically unaffected until our days, 
since the medium is transparent to their propagation. 



In addition, important improvements on the projects for direct detection of 
gravitational waves, as well as on the measurements of astrophysical phenomena af- 
fected by them have been performed or are expected for the next decade. This will 
make possible to extract constraints and to put upper limits on the magnitude of the 
generation and amplification of gravitational waves and thereafter, on the particu- 
lar model of evolution considered for the universe. Currently, the search of direct 
measurements by means of resonant bar detectors (e.g., NAUTILUS, AURIGA) has 
been followed by the development of laser interferometers (LIGO, LISA). The sensi- 
tivities to be reached for the next generation of laser interferometers detectors could 
fix these constraints on ranges still more interesting from the point of view of early 
universe theories (JlOl). Indirect constraints can be obtained from the closure re- 
lation, nucleosynthesis bound, pulsar timing measurements and isotropy of Cosmic 
Microwave Background Radiation (CMB) (|lll|). Finally, there would be imprints 



on the CMB temperature anisotropy spectrum through Sachs- Wolfe effect (|p!2|) 
and polarization fluctuations originated by scattering Thompson on a tensorially 
perturbed metric 



This study is particularly interesting from the point of view of String Theory in 
Cosmology since gravitational waves are one of the observable tracks to be extracted 
from the fundamental theory. The decoupling of gravitons from another kind of 
matter has taken place very close to the Planck Energy Scale. Thus, the shape of the 
spectrum could carry information about the physics on such energy and lenght scales. 
A number of studies computing gravitational wave production in the context of string 
cosmology scenarii has been made in the last years. ([|l^],|jl5[,[0,|]l3). Oughted 
to the intrinsic uncertainty in the cosmological model, the spectra computed in 
these works depend on free parameters coming from the earliest stages of Universe, 
far from observational knowledge. No firm free-parameter prediction is extracted, 
although the existence of a peak and end-point in the string cosmology gravitational 
wave background has been claimed (fTF 
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In this paper, we study the production of a primordial stochastic gravita- 
tional wave background generated in the context of selfconsistent string cosmology 
19| , 10 , I^TI , ) . From both treatments therein explained, low energy effective 



string action and selfconsistent solutions of General Relativity plus string matter, we 
have obtained a minimal model for the evolution of the scale factor |2^. The String 
Driven Cosmological Background is characterized by an inflationary inverse power 

2 2 

stage a/(t) ~ (tj — t) '^+^ plus a radiation dominated stage aii{t) ~ (t — tii)'^+^ and 

2 

a matter dominated stage ani{t) ~ t^. The cosmological model is linked to the Uni- 
verse observational information by means of descriptive temporal variables as per- 



formed in p4| and p3[. Our results are thus expressed in terms of standard times for 
inflation-radiation dominated transition and radiation dominated- matter domi- 
nated transition 7^. The conformal time variable rj has been defined with totally 
satisfactory continuity conditions in the transitions. 

We compute exact expressions for the power spectrum P{uj)duj and contribu- 
tion to the energy density flew of primordial gravitational waves generated in the 
transition among the inflationary and radiation dominated stages. These expres- 
sions are fully predictive and free-parameter, since the cosmological model has been 
constrainted with minimal and well stablished observational Universe information 
(the transition times %. and %n). Having fixed a minimal but satisfactory model 
for the scale factor, we study the production of gravitational waves by making three 
approaches in the perturbation treatment. 

The role of the dilaton, a characteristic field in low energy effective string theory 
and present in our model, is proved to be crucial since gives a strong signature on the 
shape and magnitude of these spectra. The computations of graviton perturbations 
with no dilaton role are showed to be equivalent to standard Quantum Field Theory 
computations, like those reported in Allen |^ and Grishchuk , but notorious differ- 
ences are found in the spectra oughted to the particular character of the inflationary 
string driven scale factor, an inverse power type. We obtain an exact expression for 
the power spectrum and energy density contribution (eqs.(^]38) and (|5.45|) ) in terms 



of Hankel functions. The power spectrum asymptotic behaviours at low and high 
frequencies are both vanishing with dependences cus and (jJ~^ respectively. This gives 
a gravitational wave contribution to the energy density asymptotically constant at 
high frequencies on Qqw ~ 10"^^. The slope change produces a peak in the power 
spectra around a characteristic frequency totally determined by the model of value 
LO^ ~ 1.48 Mhz. 

Computations with partial dilaton role, considering this in the perturbation 
equation but not on the perturbation itself, show some differences with the last 
ones, althought the general characteristics are very similar. The formal expressions 
for the power spectrum and energy density contribution are the same of the No 
Dilaton case (eqs.( ^738l) and ( |5.45D ), but they involve Hankel functions with different 



orders. Both asymptotic regimes for small and high frequencies vanish again, but 
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now with dependences cus and respectively. The peak appears around the same 
characteristic frequency, but its value is one order of magnitude lower than in the No 
Dilaton well as the asymptotic constant contribution to the energy density 

In the opposite case, when the full dilaton role is accounted, general characteris- 
tics as well as orders of magnitude of the graviton spectrum are drastically modified. 
The formal expressions are different (eqs. ( P^8^ ) , -841) ) . We extract the asymptotic 



behaviours too: for low frequencies, the dependence is vanishing again with us. 
For high frequencies, an asymptotic divergent behaviour proportional to uj is found. 
Similarly, the contribution to the critical energy density in gravitational waves is 
divergent at high frequencies, as u"^. The change of slope is less visible and no clear 
peaks are found. The transition from the low frequency to the high frequency regime 
is slower than in the previous cases (with no dilaton or partial dilaton roles) and 
full analytical expressions are needed on a wider range 10^ ~ lO^Hz. Comparaison 
of these summarized results can be found in Tables (0), (|D and (|]) and in Figures 
(i and (I). 

The behaviour obtained in the full dilaton treatment is a common characteristic 
to the String Cosmology low energy effective treatments. We conclude that the 
Brans-Dicke frame (dilaton field) involved and the type of scale factor predicted 
(inverse power type) are responsible of these features. More conclusions (the dilaton 
role, string and no-string cosmologies) are given in Section VII. 

In Section II, we recall the derivation of the background solutions obtained 
in String Cosmology, the richness of String Theory applied to cosmology is visible 
along this section. In Section III we summarize the String Driven Cosmological 
Background in order to obtain an approximative but suitable description of the cos- 
mology of our observational Universe. We translate this description to a conformal 
time-type variable more convenient for further work. Section IV shows the compu- 
tation of gravitational wave production in the framework of Quantum Field Theory 
and shows how the dilaton role can be considered. In Section V we particularize to 
our cosmological description and obtain exact expressions in suitable measurement 
units both for power spectrum and energy density in each one of the three treat- 
ments. In Section VI we analyze the asymptotic behaviours. Finally, we elaborate 
and discuss our results and compare to other String Cosmology computations as 
well as to pure General Relativity gravitational wave computations, and present our 
conclusions in Section VII. 
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2 String Driven Cosmological Background 



The String Driven Cosmological Background is a minimal model of the Universe 
evolution totally extracted from effective String Theory. Details on the model can 



be found in Ref . ||23|| , but here we summarize it. Two ways allowing extraction of 
cosmological descriptions for the background from string theory have been used. One 
is the low energy effective string equations plus string action matter. Solutions of 
these equations are an inflationary inverse power law evolution of the scale factor, as 
well as a radiation dominated behaviour. On the other hand, selfconsistent Einstein 
equations with a classical gas of strings as sources will give us again a radiation 
dominated behaviour and a matter dominated description. 

In this whole and next sections, unless opposite indication, the metric is defined 
in lenght units. Thus, the (0,0) component is always time coordinate T multiplied by 
constant c, that means t = cT . Derivatives are taken with respect to this coordinate 
t. Otherwise, constants c, % and G are explicitly showed. We consider a spatially 
flat, homogeneous and isotropic background and we write the metric in synchronous 
frame (^foo = I, goi = Q = QQa) as: 

Qf^y = diag(l, -a^(t) Sij) (2.1) 



The earliest stages in our model are provided by the low energy effective context. 
The scale factor is extracted by extremizing with respect to dilaton field, graviton 
field and matter sources the low energy string effective action (to the lowest order in 
expansion of powers of a'), which in the Brans-Dicke or string frame can be written 

as 



25|,1T|,[|T4|,|2T 



S = "Y^^ / d'^'^^le-^ [r + d,4>d^4> + (2.2) 

where Sm is the corresponding action for the strings as matter sources, H = dB is 
the antisymmetric tensor field strenght and V is related with dilaton potential and 
vanishing for some critical dimension. The dilaton field 4> depends explicitly only 
upon time coordinate. D is the total spacetime dimension, D = d+1 where d is the 
number of spatial dimensions. The string matter is included as a classical source 
which stress energy tensor in the perfect fluid approximation takes the form: 

r/ = diag(p(t),-P(t)(5,^) (2.3) 

where p is the energy density and P the pressure for the matter sources. Here, we 
will consider antisymmetric tensor H^aii as well as the potential B vanishing. We 
do not consider here the effects oughted to a non-vanishing dilaton potential. Thus, 
the low energy effective equations are obtained: 
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-2(f) + dH^ 
Hi) 





IGttGd _ 4 

— 



(2.4) 



2(H 



16nG 



D 



p 



^ and we use the shifted string duahty invariant expressions for the 
(p — \nJ\ g I, matter energy density p = paf" and pressure p = Pa'^. 



where H 
dilaton d) 



The inflationary String Driven stage appears as a new selfconsistent solution of 
the low energy effective equations ( p^ ) sustained by a gas of stretched or unstable 
string sources (|[T^,pO|). The pressure and energy density of this string behaviour 
in curved background satisfy the equation of state P = 7p as matter sources, where 
(i 7 = — 1 for the unstable behaviour in the metric here considered. Thus, the 
equation of state to be considered for the string sources is given by: 



P 



d 



P 



With this, the selfconsistent solution of effective equations (^.4]) is: 



0/ + 2d\na{t) 



ait) 

m 

pit) = pMt))^'-'^ 

Pit) ^ -I Pit) ^ 



< t < tr < ti 



-^(a(t))(-^) 



Q 



d+l 



(2.5) 



(2.6) 



Notice here t is the proper cosmic time coordinate, running on positive values; d is 
the number of expanding spatial dimensions; pi, (pj are integration constants and 
Aj, tj parameters to be fixed by the further evolution of scale factor, the parameter 
ti being greater than the exit of the inflationary stage t^. 

It is convenient for further work to express the solutions in conformal time rj 
such dr]=^: 



ttinfiv) ~ iVi - TlT 
Ad 



(2.7) 



<Pinfiri) 



ln(r77 - Tj) 



where rji is a parameter to be determined in the matching with the next cosmological 
stage. 



As discussed in ref.||23|, this is a new string cosmology solution, obtained with- 
out exploting manifestly duality relations among the solutions, and describing an 
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accelerated expansion running on the positive branch of cosmic temporal coordi- 
nate. In this way, this solution is different to the "Pre-Big Bang" solutions found in 
literature [111 . 



On the other hand, the dual to unstable string behaviour in curved backgrounds 
follows a typical radiation type equation of state P = ^ p |T^. The effective equa- 
tions (|2.4|) plus this equation of state give the following behaviour describing a 



radiation dominated stage [T^, |]l4 



An t 



R 



tr <t 



'II 



a{t) 

m 

Pit) = pu{a{t))-^'^'^ 
Pit) 



R 



d+1 



(2.8) 



i Pit) = Pf{a{t))-^'-'^ 



here 0//, 
evolution 
Again, we 
the form: 



Pii are integration constants, and An a parameter to be fixed by the 
of scale factor. Notice the dilaton remains "frozen" at constant value, 
define the corresponding conformal time to this stage. The solution takes 



dradiri) ~ (^7 - miY , r = J— J (2.9) 

(j^radiv) = 4>ll 

where rjn, (pn are constants to be determined with the scale factor evolution. 

Cosmological backgrounds can be found too, as selfconsistent solutions of the 
General Relativity Einstein equations selfsustained by the string sources evolving 
in them, as developped in |T^,[^,[^. In the curved backgrounds here considered, 
the Einstein equations take the form: 



-d{d-l)H' 



P 



(2.10) 



(rf - l)if + P + p = 



and the matter source is described by a gas of non interacting classical strings, 
whose equation of state includes the different possible behaviours of strings in curved 
spacetimes: unstable, dual to unstable and stable, each one with string densities 
V and S respectively. Taking into account the properties of each behaviour (see 
[p^ ), the energy density and the pressure of the string g described by: 



P 
P 



— ^7 iuait) + ^ + S 

{a{t))'\ a{t) J 



1 



V 



d (a(t))'' \ait) 



-Ua{t) 



(2.11) 
(2.12) 
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Equations ( [^.IID and ( ^.121 ) are qualitatively corrects for every t and become exacts 
in the asymptotic regimes, leading to obtain the radiation dominated behaviour of 
the scale factor, as well as the matter dominated behaviour. The first one is obtained 
in the limit a{t) and t 0, where the dual to unstable behaviour dominates in 
the equations and (|2.12|) and gives us: 

p{t) V {a{t))-^''^^^ (2.13) 
P{t) ^^V {a{t)y^''+^^ 

(Jj 

Solving the Einstein equations (|2.1CI| ) with sources following eqs.( p.l3| ), yields the 
scale factor solution: 

a{t) describes the evolution of a FRW radiation dominated stage, here the parameter 
tjj will be fixed by the further evolution of the scale factor. This is coherent with 
the fact that dual to unstable strings behave in a similar way to massless particles, 
i.e. radiation |TP[. In conformal time and rearranging constants, this solution reads: 



d-1 

Notice that the scale factor now described is equivalent to the before obtained in 
low energy effective treatment, eq.(p.9|). 

In the opposite limit a{t) —>■ oo, t —>■ oo, the unstable density U ^ and the 
stable behaviour S becomes dominant. The equation of state reduces to: 

p^S{a{t)y'^ (2.15) 
P = 

and from solving eqs.( |2.1(J0 with eqs.( |2.15D , the solution of a matter dominated stage 
emerges (Notice that stable strings behave as cold matter): 

here tm is fixed by the whole scale factor evolution. This scale factor in conformal 
time takes the form: 

2 

amativ) = OLHiiv - Viii)"' 1 ^ = 'JZr2 (2-17) 

here the constants am and r]jii will be determined from the transition between 
radiation and matter dominated stages. 
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The evolution model with the String Driven inflationary stage eq.( |2.6| ), followed 
by a radiation dominated stage (eqs.( p.8|) or ( |2.14]) are equivalents for our purposes) 
and a matter dominated stage eq. (|2.16|) is the String Driven Cosmological Back- 
ground that we will study. 



Notice that we consider the dilaton field remaining practically constant and 
vanishing from the exit of inflation, as suggested in the String Driven Radiation 
Dominated Solution, until the current time. Notice also the Brans-Dicke frame for 
the metric-dilaton coupling where the string action has been written eq.( |2.2| ) and 
both an inflationary and a radiation dominated stages have been extracted. The 
last one can be obtained too selfconsistently from the Einstein equations plus string 
matter, the same treatment where the matter dominated current stage is obtained 
too. In the table (|ip the scale factor, dilaton and equation of state obtained from 
these treatments can be observed in a comparative way. 





String Driven Solutions 


String Sources in G.R. 




LEE + {U,V) strings 


EE + gas (U,V,S) strings 




inflation 


radiation 


radiation 


matter 


a{t) ~ 


(ti - 1)~^ 


2 

td+1 


2 

td+1 


2 
t-d 


a{r]) ~ 


2 


2 


2 

iv - r]ii)d-i 


2 

(?? - mu) ""-^ 






cte 






0(r/) ~ 


(-^) Hvi - V) 


cte 






p(t)~ 


(a(t))(^-^) 




V{a{t)y^^^'^^ 


S{a{t))-' 


string 


unstable 


dual to u. 


dual to u. 


stable 


sources 


P = -\P 




P=\P 


P = 



Table 1: Selfconsistent solutions giving the String Driven Cosmological Background 

Table comparative of the solutions for the scale factor a{t), a{rj). t is cosmic time, r] 
is conformal time and (j) is the dilaton field. LEE means low energy string effective 
equations in the Brans-Dicke frame, for a (d-f- 1) homogeneus, isotropic and spatially flat 
background; antisymmetric tensors and dilaton potentials are neglected. EE means the 
General Relativity Einstein Equations plus a gas of strings as classical sources matter. 
The dominant string behaviour as matter sources and their corresponding equation of 
state are written in the last row. 



3 Scale Factor Description 



Since the behaviours above extracted are asymptotic results, it is not possible 
to give here the detail of transitions among the different stages. It would be ex- 
pected some changes in the regimes or equations governing each stage and leading 
to the next one, but this is an open question in the framework of string theory both 
for inflation-radiation dominated as well as radiation dominated-matter dominated 
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transition. Taking the simplest option, we consider the "real" scale factor evolution 
minimally described as: 

ai{t) = Ai{ti - ty'^ teiU,tr) (3.1) 
aii{t) = Ajit^ te{tr,t^) 
ajjj{t) = Am t"" teit^,to) 

with transitions at least not excessively long at beginning of radiation dominated 
stage and tm beginning of matter dominated stage. We define also a beginning of 
inflation at ti and a current time to- 



From the point of view of String Theory applied to Cosmology, it would be 
reasonable to have not instantaneous and continuous transitions at tr and l|2B 



But this condition greatly difficults the computation that we want make because 
it introduces free parameters, loss of physical meaning and effects of an unknown 
transition dynamics 0] oughted to discontinuity on the scale factor. In order to 
construct a suitable minimal model for gravitational wave computations, we will 
merge our lack of knowledge on real transitions by means of descriptive temporal 
variables in function of which the modelized transitions at ti and ^2 are instantaneous 
and continuous. 



We link this descriptive scale factor with the minimal information about the 
observational Universe. We consider the standard values for cosmological times: the 
radiation-matter transition held about 7^ ~ lO^^s, the beginning of radiation stage 
at Tr ^ 10~^^s and the current age of the Universe % ~ Ho~ ~ 10^'^s (The exact 
numerical value of % turns out not crucial here). We impose also to our description 
to satisfy the same scale factor expansion (or scale factor ratii) reached in each one 
of the three stages considered in the real model (pTT|). Explicit computations can 
be found in |^ and leds finally to this scale factor written in cosmic time-type 



variables in the more convenient way for our further work: 



di{i) = Aj{ti -t) ^ U<t<ti (3.2) 
an it) = An {t-tnf h <t< h 



ani{t) = Ajnitf t2<t< to 

with continuous transitions at ti and ^2 for both the scale factor and first derivatives 
with respect to the descriptive cosmic time-type variables i, i and t. In terms of the 
standard observational times tr and tm, the transitions ti, t2 and the beginning of 
the inflationary stage description ti are expressed as: 



M V MJ 
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and the set of parameters of the scale factor ( |3.2| ) can be written too in terms of tr, 
tm and the global scale factor An: 



Notice that the intermediate descriptive scale factor could not be physical, but 
this is not important under the purposes here. Since the description is linked with ob- 
servational information and it satisfies the proper scale factor ratii, our computation 
will be equivalent to that made on a full and physical model. These are equivalents 
unless the features left on gravitational wave production by the details of an in any 
case unknown dynamics of transition and minimized through their modelization in 
satisfactory continuos way, as made here. The time variables t of inflationary stage, 
and t of radiation stage are not a priori exactly equal to the physical time coordinate 
at rest frame (multiplied by c), but transformations (dilatation plus translation) of 
it. The low energy effective action equations from where the scale factor, dilaton 
and equation of state have been extracted, allows this transformations. With this 
treatment of cosmological scale factor, we will attain computations free of "by hand" 
parameters, and with full predictibility as can be seen in the next sections. 



It is also useful to express this model in conformal time variables 77, defined 
as drf = We will have a conformal time variable for each stage (fj, fj and 77) 
oughted to the different scale factor shape and different descriptive temporal variable 
on which we integrate. In this way, the conformal time variable constructed enjoys 
continuity on their derivatives along the transitions at ti and t2 [24]. 



dr] 

dt 



dr] 



dr] 



t2~ 



dr] 
'dt 



t2^ 



We can impose also continuity conditions in order to have a conformal time uniquely 
defined at transitions. r]i = f]{ti) = f]{ti) and r]2 = 77(^2) = ^(^2) (In this respect. 



different ways can be found in literature, see We use the remaining freedom in 

put a simplicity condition, asking a simple conformal time dependence in radiation 
dominated stage. It yields the next expressions for the scale factor in conformal 
time: 



aniv) = anivY 
aiiiiv) = o^iiiiViii + vY 



r]i<f] <r]i 
r]i <f] <r]2 
r]2<r] 



(3.5) 



with continuous and suddenly transitions (continuity of scale factor and first deriva- 
tive with respect to conformal time) at r]i and 772- The six parameters satisfy these 
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matching relations: 



\r+q 



Vi ^ Vi U- + - 



Oiiii = ajj(^y^ {r]2) 



(3.6) 



Because its construction on a cosmic time description with fully continuity condi- 
tions, all the parameters involved in this conformal time scale factor have a complete 
and unique expression as functions of the observational standard transition times 
tr and tm, beginning of inflationary stage tj, the exponents Q, R,M of cosmic time 
dependences and the global scale factor constant An: 



Q 



R 



M 



Q + 1 



R + Q 



1 - R 



m 



1-M 



R(Q + 1){1-R) An 



R' 




^ 1 


M. 




{Q + if 


"(1 


- R)^ 





R+Q 



Q+l 



M-R 



Vni = 



R(l-M)(l-R) An 



am = An{l - M) 



M 



M 



R-M 



(3.7) 



Meanwhile, the conformal time transitions rji, 772 and rji have the expressions: 

l-R / „ N l-R 



Vi 



(1 - R)An 

l-R 



{Q + l)An 



(1 - R)An 

R + Q _ Q f U-ti 
R{l-R) R\tr-ti 



(3.8) 



It is also possible give expressions for the conformal time in each stage as function 
of the proper cosmic time: 



V 



l-R r 



An{Q + l) 



R + Q Q( ti-t ^'^+^ 
R{1 - R) R \ti - tr 



An{l - R) 
An{l - M) 



R-M / t ^ 



R{1 - R) R \tr, 



i ^ (^2 5 ^y) 

t G (tri ^m) 

t e {tm, to) 



(3.9) 
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Finally, it is useful for further work to compute here the value for the conformal 
time and scale factor at current time r]o. If we define to = cTq and from eqs.(|3^), it 
is easily seen: 



^0 



,1-R 



Ajj{l-M) 



R-M M /to^ 



+ 



R{1 -R) R \t 



(3.10) 



The scale factor written in conformal time for current time takes the form a///(?7o) = 
OiiiiiViii + Vo)"^- From eqs. (|377|) and ( p.lO| ) we translate this value to observational 
transition values: 



aiiiiVo 



R 



4. M J. R-M 



(3.11) 



The last point is to make an approach for the dilaton field. This is considered 
practically constant from the beginning of radiation dominated stage until the cur- 
rent time. Their value can be supposed coincident with the value at exit inflation 
time in a sudden but not continuous transition for the dilaton, since no one of their 
temporal derivatives can match this asymptotic behaviours. 

hi = <Pivi) = 01 (3.12) 
Remember the expression for dilaton in inflation dominated stage, that gives 

= (Pi + 2dlna{r]i) (3.13) 



temporal dependence 
parameters 


Strir 
d 


ig Driven Cosm.B. 
d = 3 


Cosmic 
Time 
description 


infl. Q 


'2 
d+1 


1 

2 


rad. R 


d+1 


i 

2 


mat. M 


2 

d 


2 

! 


Conformal 

Time 
description 


infl. q 


•2 

d+3 


3 


rad. r 


2 

d-1 


1 


mat. m 


2 

d-2 


2 



Table 2: Temporal Dependences in the String Driven Cosmo logical Backgrounds 

Values of the exponents in the temporal dependences for the inflationary, radiation dom- 
inated and matter dominated stages in String Driven Cosmological Background. The 
parameters Q, R and M appear in the cosmic time description, whereas the parameters 
q, r and m are in the conformal time one. Particularization to the three dimensional 
case is given. 
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4 Generation of Gravitational Wave Perturbations 
in String Cosmology 



After obtaining the minimal cosmological model for the scale factor in the 
framework of string cosmology, we will compute one of its important observational 
consequences, the stochastic background of relic gravitational waves. As first noticed 
by Grishchuk 0, relic gravitational waves are inevitably generated in practically all 
cosmological models oughted to their variable gravitational field. 



The generation of gravitational wave backgrounds can be studied classically 
as propagation and amplification by the expanding background of one initial spec- 
trum of tensorial perturbations on the metric. From the quantum mechanical point 
of view, a particle production mechanism occurs and gravitational waves of every 
wavelenght are generated as gravitons from an initial vacuum state. This translates 
into a combined classical-quantum mechanical formalism. (See for instance 0] and 



2^). Gravitons produced at scales near the Planck time must exist today in the 
form of gravitational waves of wavelenghts to be determined by the propagation 
equation. 



The number or particles produced between an initial -asymptotic- vacuum state 
and a final state, where particles -gravitons- has been created by the effect of the 
changing gravitational field can be obtained by computing the Bogoliubov coeffi- 
cients 0]. Grishchuk has a similar approach modelizing the changing gravita- 
tional field as a potential barrier. Each polarization degree of freedom of tensorial 
perturbation satisfies the same propagation equation that a massless scalar field 
(0,|Q,0). But in frameworks with a dilaton field, attention must be put in ac- 
counting for the dilaton contribution to the gravitational field change and to the 



metric perturbation, breaking the validity of reported treatments |15 



In the next, we summarize and follow the treatment of tensorial perturbations 
in cosmological backgrounds within the Brans-Dicke frame developped in ref. [l^ . 
That is necessary for accounting of the full effect of the dilaton field and for the 
consistency with the framework in which the string cosmology solutions have been 
obtained. We will show that this formalism includes as a particular case (no dilaton 
field) the formalism developped by Allen ^ and Grishchuk in the framework of 
quantum field theory. 



The general action written in Brans-Dicke frame takes the form: 

3 

S = -t/t^ / d^^'x^\e~^ (R-uj g^^ d^<P d,cj)) + Su (4.1) 
IottG/) J ^ 

where Sm represents the contribution of matter sources such that its stress tensor 
IS J\ g \T^i, = 2|^, antisymmetric tensors and dilaton potential are vanishing and 
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uj is the usual Brans-Dicke parameter. For = oo General Relativity expressions 
are found and for uj = —\ this action coincides with the low energy effective string 



action eg . ( ^I2|) . By varying with respect to the metric g^^y and to dilaton the 



following two equations are obtained: 

i?/ + V^V^^ + + 1)[(5/((V0)' - □(/)) - V^0V>] = ^^e^T/ (4.2) 

+ cu(V0)^ - □(/) = (4.3) 

where □ = f^'^V^Vi^ is the covariant D'Alambertian operator. Again, the case 
= — 1 recoveres the low energy string effective equations ( |2.4|) . 

The free- linearized wave equation for a generic metric fluctuation h^^^y = 5g^v is 
obtained by expressing the perturbed metric as g^u = giiv + ^giiu, and all the sources 
are considered fixed 5T^" = 50 = 0, second order (and higher) terms in h^y are 
neglected. 

In the following, the study is made only for the tensorial component of the 
metric fluctuation h^^. The general case for the metric is considered: an anisotropic 
background, with d expanding external spatial dimensions with scale factor a{t) and 
n contracting internal dimensions with scale factor 6(t), both fiat and Euclidean. 
Expressed in the synchronous gauge, the metric is: 

g^, = (1, -a2(f)7,,(x), -h\t)-iah{y)) (4.4) 

where again the component (0,0) is written in lenght units and corresponds to the 
variable t = cT, indexes i, j run from {1 . . .d) and indexes a, b run from {d + 1 . . .d + 
n+l). 

The metric perturbation is considered propagating only in the external spatial 
dimensions, decoupled from the sources and purely tensorial: 

h^^ = h^^{yi,t) (4.5) 

V = V = (4-6) 

That means, a pure gravitational wave that can be expressed in the transverse- 
traceless gauge, thus g^'^ h^^ = and V^h^" = 0. With these conditions, eq. (|4.3| ) is 



trivially satisfied and from ( ^.2| ) the free-linearized wave equation for the tensorial 
metric perturbation h^^i, is obtained: 

+ i V g"'' - /i''"V^V„0 = (4.7) 

It must be noticed that this equation is cu-independent and therefore it is valid in 
general for all Brans-Dicke cases and in particular for the low energy effective string 
action. Eq. ( [4.7| ) takes a shorter form when the whole information about the metric 
is introduced [p!5| . 

nhi^ - (Phi = (4.8) 
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Writting the perturbation in terms of eigenstates of the Laplace operator of 
wavenumber k: V^hi\k) = —k'^hi^{k), eq. (|4.8| ) takes the form: 



hJ + {dH + nF- (f))hi' + - V = (4.9) 



a . 



here dot means derivative with respect to t (remember t = cT), H = ^ and F = |. 
The most suitable form of eq. ([4.9|) is obtained by writting it in conformal time 



df] = -%r. It is convenient to define a variable if): 



ijji^ = hi^ 6^ e"2 , (4.10) 



in such a way that is identified with each one of the two polarization modes of 
the perturbation hij. From eq.( |4.9| ) it is immediately seen that each polarization 
mode ifj satisfies the corresponding equation written in conformal time: 

^/j" + {k"^ - V{r]))ij = (4.11) 

where : 

1/2 1 , b 1,, n (b\ I 



Thus, the problem of metric fluctuation propagation is reduced to a second 
order differential Schrodinger type equation, where the "potential" V{ri) takes into 
account the conformal time derivatives of the background as well as those of dilaton 
fields; that means, the whole variable gravitational field "pumping" the generation 



and amplification of gravitons.[15 



Following our definition of the variable t = cT in length units, the confor- 
mal time variable t] has also length units. Thus, the "potential" ( [4.12| ) have units 
(length) and the wavenumber k has units (length) 



As it was pointed out in ]T5|, this treatment generalizes to higher dimensional 



backgrounds the treatment more widely work, usually based on the pioneering works 
of Grishchuck [ffl], pi. It generalizes also other early works like those of Allen M. 
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4.1 Particular Cases 



String Driven Case 

In order to compute perturbations in our cosmological background, we will 
consider the suitable particularization of eqs.( [4.11| ) and ( ^.121 ). For arbitrary d and 
77, = eqs.( [4.1l| ) and ( 4.12|) take the form: 



ip + (r - V{r]))ilj = (4.13) 

+ ,4.14) 

The perturbation variable ip is written as: 

tpi^ = hia^e~^ , (4.15) 
and as consequence, "00 = V^)! = 0- Equivalently, the tensorial perturbation is: 

(4.16) 





= ipi^a 


d+3 4> 


2 62 


ho fi 


= h^Q 


= 


ha fi 


h^ a 


= 



where the index /z runs now from 1 . . .d. 
No Dilaton Case 



It is possible to extract from eqs. (|4.13|) , ( [4.14| ) and ( [4.15| ) the case without 



dilaton field, in which we obtain the set of equations: 

/ + (fc^ - V{r]))ij = (4.17) 

The perturbation variable i/j is: 

^/ = hi a'^ (4.19) 

and the metric tensorial perturbation is: 

hij = i/ji^ 

ho fi = h^Q = 
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The case of three spatial dimensions gives the simplest expressions: 

i)" + - V{7]))i) = (4.20) 
V{V) = (4.21) 

where the perturbation is simply written 

iPi' = a (4.22) 

with hij = ipi^ and ho ^ = h^Q = 0. This particular (3+l)-dimensional 
case gives the same set of equations ( ([4.20|) - (|4.22|) ) that can be found in Grischuhck 
treatments [Q. These equations can be also identified with the formalism found in 



For the sake of completeness, we compare now with the linearized Einstein 
equations for an homogeneus, isotropic and spatially flat metric gf^^, which in the 
conformal time gauge is: 

ds"^ = a{r]f{-dTf + dx^) (4.23) 

A gravitational perturbation h^^ = a'^{r])f]iu of comoving wave number k is 
overimposed such that: 

9,.u = a\r]){r]^, + f^,) (4.24) 



By performing a Fourier expansion, f^^, can be described as: 

f^^uiv,^) = ¥'(??) e/,i.(k)exp(zkx) 



(4.25) 



where e^j^ is the constant polarization vector. It can be shown that the amplitude 
perturbation satisfies the equation M: 



a{r]) 



(4.26) 



where =| k |, k is the comoving wave number and prime stands for conformal 
temporal derivatives. This equation is equivalent to eq. ( [4.9| ) in the Brans-Dicke 
formalism, provided particularization to isotropic cases {F = 0) and no dilaton 
dynamics {(/) = 0). For d = 3, eq.( [4.26D is the known equation found in the literature 



Let us write: 
Then, we have 

u"{r]) 



f d-i \ 



d-l(3-d(a'i7^)y a"iv) . 

2 1, a{v) J a{v) ' ^ 



iL{r]) = 



(4.27) 
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that, in the particular three spatial dimensions take the form: 



a{r]) 
a{7]) 



u"{r]) + - V{7]))u{r]) = V{r]) - """^"^^ 



a{r]) 



This shows the equivalence with the Grishchuk formalism and with the (3+1)- 
dimensional particularization of the No Dilaton Case showed above. 



4.2 The Power Spectrum 



We perform an expansion of the vacuum state (amplitude of perturbations at 
the classical level) in positive- and negative- frequency modes in each stage. The 
Bogoliubov coefficients mix the modes before and after the transition or the graviton 
creation- annihilation operators and ak'^. \ f3k l"^ gives the number of particles of 
comoving frequency k generated on the vacuum state. This is the "initial" state 
(at a time sufficiently early before the transition to radiation dominated stage), so 
the corresponding Bogoliubov coefficient in the radiation dominated stage yields the 
total number of particles created during the whole inflationary era with respect that 
initial state. 



We will determine the Bogoliubov coefficients a{k) and (3{k) for every wavenum- 
ber perturbation k by considering a continous matching of the perturbation ampli- 
tude at the transition time rji between inflation and radiation dominated stages. 
Thus, we will impose continuity of the amplitude perturbation -treated as quantum 
state- as well as of its flrst conformal time derivative at the sudden transition at 
conformal time r]i. This is similar to work with the picture of a potential barrier 
generated by the variable gravitational fleld 0. Gravitational waves would inter- 
act with this barrier. By describing one state "out" the barrier, one state "in" 
with frozen amplitude and a comeback to "out" , the Bogoliubov coefficients can be 
computed. But as above said, the contribution of the dilaton change can be missed. 

The stochastic background of gravitational waves generated in the cosmological 
transitions is described by the power spectrum P{uj)duj, which is the energy density 
contained in gravitational waves with frequencies in the interval between u and 
u) + du; oj is the proper frequency at the detection time related to the comoving 
wavenumber k through: 

a. = ^ (4.28) 
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The total energy density in gravitational waves is computed as pcw = / P{uj)du 
[|]. The power spectrum P{u)du is given by P{uj)du = 2hu!\ j3 fdudN where dN is 
the density of states, which in three spatial dimensions (as we have at the detection 
moment) is dN = ^^^du; \ P f is the number of gravitons created in each frequency 
interval and each polarization state. This comes multiplied by the two polarization 
degrees of freedom of the tensorial metric perturbation and by the energy hcu of each 
one: ^ 

P{u)diu = ^\pfduj (4.29) 



Thus, the power spectrum P{uj) has dimensions [-P(u;)] 



energy time 



volume time lenght ' 



It is usual to write the power spectrum as fraction of the critical energy density 
by octave, which is the closest form to gravitational wave detection procedures: 

_ 1 dpGw{uj) 

I Low 



Pc dlnu 

here pcw is the energy density contained in gravitational waves and pc is the critical 

SvrG 



energy density pc = This is equivalent to compute 



now = - P{^) (4.30) 

pc 



Obviously, VLqw is dimensionless [^gvk] = 1- 



5 Power Spectrum in String Driven Cosmology 



5.1 No Dilaton Case 



We begin our computations of power spectrum of gravitational waves with the 
simplest case: we will neglect in this the whole dilaton role. Thus, we will solve 
eqs. (|4.17|) , (^4.18|) and the reduced amplitud perturbation given by ([4.19|) . 



by: 



In the inflationary stage, from eqs.(|3.5|) and ( |4.18| ), the "potential" V{r]) is given 



d-1 



2 -i^^^-^)" 



(5.1) 



The general solution of eq. ( [4.17|) with the potential (|5.1|) is given in terms of Bessel 
functions: 

^ = (^/ - v)^{c,n,^'\k{vi - n)) + C2nJ'\k{rii - fi))) (5.2) 
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with index u: 

z/ = ±i((d-l)g+l) (5.3) 



(5.4) 



Or, from eq.( |2.7|) in terms of the number d of the spatial dimensions. 

3d+l 
~ 2{d + 3) 

It gives = I for three spatial dimensions. 

The coefficients Ci and C2 are fixed by the choice of initial or boundary condi- 
tions, i.e. by the inflationary vacuum state Thus, the asymptotic expression of 
( |5.2|) for early times should contain only positive frequency modes. 



The argument of the Hankel functions involved can be expressed in terms of 
the cosmic time. From eq.( 



z = k{r]i-i]) = k (^r]i{l + ^) - r]j , (5.5) 

and with the expressions relating the conformal time description with the observa- 
tional transition times and cosmic time eqs.(P?7|), ( |3.8| ) and (|3.9|), we get: 



k — L 



R Q 



M'J R{Q + l)An 



ti - 1 

tj — tr 



Q+1 



(5.6) 



As discussed in section II (see eq. ( p.6|) and comments below), the parameter tj is 
always greater than the exit of inflationary stage tr, thus the argument z is real and 
2; > for all t < tr- For time very early in the inflationary stage t —>■ ti and t << tr, 
the argument z approaches its maximum value. Thus, we must identify the large 
argument regime with the vacuum state at enough early inflationary stage (usually 
refered as well below the transition to radiation dominated epoch). 

The asymptotic behaviours of Hankel functions [27|: and the positive frequency 
mode condition at early stages, flx C2 = 0. The remaining freedom in the coefficient 
Ci is linked to the normalization of the vacuum inflationary state and we include it 
in the normalization coefficient A/". In this way, the reduced amplitude perturbation 
in the inflationary stage (before the transition) is written as: 

i^i = ATivi - vY^'H^'^ iKvi - m (5.7) 



In the radiation dominated stage, from the scale factor eq.( |3.5| ), the potential 



( 4.18 ) takes the form: 
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And the general solution of eq.( [4.17D , with eq.( |5.8D is given again in terms of 
Hankel functions: 



but now with index /i: 



I, = ±^iid - l)r - 1) . 



Since the radiation dominated stage has r = (eq.( p^.9| )), it makes: 

1 

/^ = ± o 



(5.9) 
(5.10) 

(5.11) 



In this case, by using the properties of Hankel functions with fractional index (See 

ion: 



, ijj eq. (|5.9|) has the expression: 

i/j = —I 

which can be written also in the way 



V TTK ^ 



(5.12) 



(5.13) 



Here we have defined Ai a suitable normalization constant and the coefficients 
a = ^ g- Q^Yid p = g«fcm are the Bogoliubov's coefficients, since they are 
multiplying the positive- and negative-frequency pure modes. 

With the expressions for the reduced amplitude perturbation in each stage ( |5.7| ) 
and ( [5.13[ ), where we remember the parameter u has the value given by (|5.3| ), we will 
determine the value of Bogoliubov's coefficients a and (3 by matching the reduced 
amplitude perturbation and its first conformal time derivative continuously across 
the inflation-radiation dominated transition value 771. 



drj 



di/j 



II 



V = Vl 



drj 



(5.14) 
(5.15) 



V = Vl 



The derivative of the reduced amplitude perturbation can be computed as: 



djJi 
dfj 



— Mk{'r]i — rf) 



d{k{rii - f])) 



^ = iM\ \({-ik)ae-"'^'^-''''^ + (iA;)/?e''=(^-''^) 
di] \ Tck^ 



(5.16) 
(5.17) 
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In the next, we will call: 



(O 



d{{k{r]i - 1]))) 



(5.18) 



7y=r)i 



In particular, the argument of the Hankel functions involved in the matching 



can be expressed, following eq. (|3.6|) 



Kvi - v) I 



V = Vi 



k7]i^ = X 



(5.19) 



With this property, equations (|5.14|) and ( ^.15|) can be written as: 

Ar(X)^7^,(i) (X) = iM^ia + /3) (5.20) 

^Ar(x)" 5 7^:^(1) (X) + kUixyn'u^^^ {x) = mV2^{^ - a) (5.21) 

2k 



By solving this system, we find for the Bogoliubov's coefficients: 



2 2A^' 
^ V 2 2A^ 



\2X 



(5.22) 
(5.23) 



From the properties of the Hankel functions (See [^), and their conjugated, 
for real p and X which is the case here: 



n,^^\x) = n,^^\x) 
T(jJ\x)=n',^^\x) 

The conjugate of Bogoliubov's coefficients a and J3 can be written as: 



^ A/" , 1 

OL= — \l— -r^X'^ 



2 2M' 



i^_,\n,m{x \ + -Hj^Hx 



\2X 
1 



(5.24) 
(5.25) 



(5.26) 
(5.27) 



With the properties of Hankel functions (See |]28|) and their wronskian (See 
): 

4 



nz 



(5.28) 
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and by rearranging coefficients, we obtain the expressions for the square moduhi of 
Bogohubov's coefficients. 



,2 'n\MHy-\ 



I 



V (-1 

2 + 



n,^^\x)n,^^\x) + 



2X 



j^n.^'\x)n.^,^''^{x)-i 



V 



71 [ly - -KiU 



AX 



|2 TT 



2 r 



+ 



X^ 



1 + 



X' 



n,^'\x)n,^-'\x) + 



V 



(2)/ 



AX 



u — 



(5.29) 



(5.30) 



By requiring the Bogohubov's coefficients satisfy the usual normahzation con- 
dition: 

we have the constraint: 



M 



(5.32) 



This condition holds for all fc, thus it ensure us the absence of isolated depen- 
dence of k through the normalization constants. The difference among them is only 
a phase. The same result can be reached by computing explicitly the constants M 
and M. in the framework of a specific normalization scheme. 



5.1.1 The Power Spectrum in Physical Units 



Now we are able to compute the power spectrum of the stochastic gravitational 
wave background with full physical meaning. In terms of the proper frequency io at 
detection time eq.( [4.28D , we read the variable X eq.( |5.19| ) as: 



X 



.rj c 



(5.33) 
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Notice that the dependence with the global scale factor parameter Ajj (the only 
quantity unknown in these expressions), is cancelled in X by the product of r)ia{r]o). 
With the help of eqs. (|3.8|) , (|3.11| ), ( p.7|) and some rearrangements we have: 

X = ujS , S = —-^—-(^yU^YTr (5.34) 



{Q + 1)M\T^J \% 

That is, we find the argument X related to the physical proper frequency in 
a way totally determined by the cosmological model studied, and linked to obser- 
vational times. The ratii and are exactly the scale factor expansion 
reached during the radiation dominated and matter dominated epochs, respectively, 
introduced as part of our linking to observational Universe information on the min- 
imal cosmological model. (See section III). Both values are strongly constrained by 
the observational evolution data. Meanwhile, the value % (beginning of radiation 
dominated stage) gives us the order of magnitude of the exit of inflationary stage, 
related to the order of magnitude of GUT scales. 

Obviously, the quantity S has dimensions [S] = time. Thus, as [u] = rad (time)~^, 
X is dimensionless. We can give here the value of S for our cosmological scale fac- 
tor. From String Theory, we have obtained the scale factors with parameters of 
temporal dependence Q = R = and M = | for inflationary, radiation 
dominated and matter dominated stages, respectively (See eqs.( |2.6| ), (p.8|), ( |2.16| ) 
and tabled^)). By linking with the cosmological observational description, we have 
considered the standard values for 7^, %n and Tq given in previous section. This 
leads to the expression 

Xstr.n, - U ^- j ^- j % . 

For three spatial dimensions, and the time constants given in seconds, we have: 

Xstring = ^ "10 3 seg (5.35) 

If we define a characteristic frequency which makes unity the argument X 
S = 1), the value of it would be 

uj^ ~ 9.28 10*' — ~ 1.48 Mhz (5.36) 

s 

Now, we compute the full and exact expression of the power spectrum for the 
stochastic background of gravitational waves produced in the inflation-radiation 
dominated transition. Following the usual definition, the power spectrum of gravi- 
tational waves computed in three spatial dimensions is: 

P{uj)du = -^iu^dul (3 1^ (5.37) 
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From eq.( |5.3(J| ) and eq.( |5.34| ), we have: 



P{uj)duj 



S 



-duj 



(coSf 



,(2)f 



2ujS 



where S is given by eq.( ^.34| ). 



TT V 



TT 



(5.3^ 



The expression in square brackets is dimensionless. All physical units are found 
in the term multiplying it and proportional to u'^duj. For our scale factor in three 
spatial dimensions, the coefficient in front will take the value: 



1.603 10 



-54 erg seg'^ 



Sttc^ S cm-" 
thus the power spectrum P{uj)duj has dimensions of energy density. 



(5.39) 



Notice the dependence of P{oj) on the scale factor parameters through the 
parameter i/, eq. (|5.3|) , characterizing the expression of the power spectrum and its 
uo dependence. This is the way the cosmo logical theory considered is imprinted on 
the gravitational wave spectrum. The parameter X eq.( ^.34D , is uniquely fixed by 
the scale factor and according with observational values. We will see in the next 
sections as the different dilaton roles and the different theories, although involving 
the same scale factor, led to different power spectra shapes. 

It is also possible to express the power spectrum in terms of the Bessel functions 
JAz) and the Neumann functions yuiz) 



with the Wronskian[27 



{Juiz), y^z)} = — . 

TTZ 



(5.40) 
(5.41) 

(5.42) 
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The power spectrum is expressed as: 

, 2 

-duj 



P{uj)diu = - — -uj 



n 2^ 2 



S 



^-2 



+ 



/"''^^^ (ju-iiujS)x^,{iuS) + y,^,{iuS)y,^,{LuS)) + 



2uS 



2 

AujS 



/ \ 

'^I^MuS)x^iiujS) + y,iujS)y,-,iLuS)j + 



(5.43) 



Another expression for the power spectrum is in terms of only Bessel functions 
Jv, J-v In order to make this, we use the relations (||28[, [p7[| ): 



sin(7rz/) 



smi^TTZ/ 
2 sin(7rz^ 



We find the next expression dealing only Bessel functions: 
P{u)duo = - — -cu^ „ . , duj 



1 + 



8nc^ S sin^fyrz/ 



- 2 cos{nv) J^{L0S)J-y{uS) j + 



+ 2cos(7rz/)J'^_i(a;5)Ji_^(c<j5) ^ + 



2c^5 
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2 uS cos (iriy) 



4 uS sin^ fyr i^) 



TT 



(5.44) 



Eqs. (|5.44| ) and (|5.43|) manifest explicitely the real character of the mathematical 
expressions involved in the power spectrum. In fact, the imaginary terms belonging 
to Hankel functions cancel exactly among them and with the free imaginary term. 
This is more easily seen looking in the last expression, where the power spectrum 
is fully expressed in Bessel functions (real functions when the argument is real) 
(See ||2^) From the numerical point of view, expressions (|5.43| ) or (|5.44| ) are more 
convenient, since standard algorithms are limited in compute Hankel functions to 
very high precision at large arguments and it leds to wrong results in computing 
the above expressions. Expressions ( ^.43 ) and ( 5.44 ) allow extend 1-2 orders of 
magnitude more the numerical computations of the exact spectrum. In any case, in 
the next chapter we will overcome this limitation with analytical expressions for the 
regime of small and large arguments. 



Finally, we compute the fraction of critical energy density in gravitational waves 
by octave. Following the definition eq. ([4.30|) and using the expression (|5.38| ), we 
have: 



KG [T^ 



GW 



+ 



3i/o'c5 S 



1 + 



{uSf 



^^7^._i«(^5)7^._i(^)(^5) + 



^-2 



n,^^\ujS)nJ^\ujS) + 



TT U 



TTl/y- 2 



(5.45) 



Again, we will compute the value of the coefficient multiplying u"^ and the 
expression in square brackets: 



hG 
3^0^ 



S 



2.488 10"^^ seg" 



(5.46) 
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5.2 Partial Dilaton Case 



In this case, we will consider the role of the dilaton field in the equations for 
the reduced amplitude perturbation tf). We will see that this procedure still means a 
partial consideration of the dilaton role on the gravitational wave spectra. By now, 
we will solve eqs.( [4.13D and ( [4.14|) with the reduced amplitude perturbation given 



by (EH). 



Now, in the infiationary stage, we take into account not only the scale factor 
but also the dilaton field. Following eqs.(|2.6|) and ( |3.5| ) we have for the dilaton in 
this stage: 

,j){ff) = (l)i + 2d\Yia{r}) (5.47) 



Thus, the potential eq. ( |4.14| ) takes now the form: 



'd+l 1 



Then, we have again an equation whose solution is expressed in terms of Bessel 
functions, but with the parameter v having the value: 

(d+l 1\ 
. = ±^— (5.49) 

Or, expressed in terms of the spatial number dimensions, 

d-l , , 

V = — 5.50 

2{d + 3) ^ ' 

that is, = I for three spatial dimensions, and fixing the positive sign in the 
eq.(0|). 



Notice the difference with the parameter v eq.( |5.4| ) obtained in the No Dilaton 
Case. The difference among both cases is: 

^Ipart.dil. ~ '^Inodil. + ~ 

It leds to a totally similar expression to eq. ( p.7| ) for the solution of the reduced 
amplitude perturbation: 

^/ = Mim - W'H,'^^^ ik{rji - ^)) (5.51) 
where u is given now by eq.( ^.49D and A/" is a normalization constant. 
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For the radiation dominated stage, the scale factor is given by eq.( |3.5| ) and the 
dilaton field, according to eqs.(p78|), remains "frozen" in its constant value at the 
exit inflation time: 

hi = m)\-, = ,, (5.52) 



Thus, we recover exactly the same expressions for the potential eg. ( |5.8|) , the 
general solution eq. (|0| ), the parameter /i = ±| and consequently, the same expres- 
sion in terms of Bogoliubov's coefficients a and (3: 



V tt/c ^ 



ik{r]-rii) 



(5.53) 



with Ai a normalization constant. 



The Bogoliubov's coefficients computed now have exactly the same formal ex- 
pressions yet obtained in the last section, eqs.( |5.29| ) and ( |5.30| ), since the reduced 



amplitude perturbations involved in the matching at the transition at rji are for- 
mally the same functions ( p.T] ) and ( p.l3| ) in both cases. Only the parameter i/ is 
different, now given by eq. (|5.49|) . 



Similarly, all the relations in order to compute the power spectrum eq.( p.3(j| ), the 
fraction of critical energy density flow eq.( 5.45|) and the expression of the argument 



X eq.( |5.34] ) as function of the proper frequency u, hold in this case. Notice that the 
argument X does not depend on the dilaton role; this quantity comes fixed from 
the evolution of the scale factor. Thus, we have again P{uj)duj and VLqw as given by 



eqs.( |5.38D and ( [5.45| ) but with the parameter p given in this case by eq.( p.49D 



In three spatial dimensions, we have the same values for the parameter S 
eq.( |5.35| ) and for the characteristic frequency Ux eq.( ^.36| ). The numerical values of 



coefficients multiplying the square brackets expressions in ( |5.38|) and ( p.45| ) would 



be, in principle, differents since they depend of the parameter f. But in our case, 
in three spatial dimensions, the value of (z/ — |)^ is the same for both cases, giving 
the same numerical value for these coefficients too. 

It must be observed that the normalization constants M and M. are constrainted 
by the normalization condition of Bogoliubov's coefficients to satisfy the same re- 
lation found in the case before: N and differ only in a phase. (See eqs. (|5.32| )). 



But the case of particular and explicit normalization applicable to the no dilaton 
case is not generalizable to the dilaton case. Now, both constants will remain as 
undetermined parameters in the amplitude perturbations. 
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5.3 The Full Dilaton Role 



In this section, we will take into account the full dilaton role. In the last section, 
we have considered it through the reduced amplitude perturbation ^, solution of 
eq.( |4.13| ). But we have left unconsidered the dilaton role in the total amplitude 
perturbation /i^j,. In order to include this, we will match in this case the metric 
amplitude perturbation h^y instead of only the reduced amplitude perturbation tfj. 

In the framework of General Relativity, as the treatment of Grishchuk |^ and 
Allen [|], both procedures are equivalents, since the metric perturbation and the 
reduced amplitude differ only through a power of the scale factor. In those cases, and 
with a description with suitable properties of continuity for the scale factor and its 
first conformal time derivative like those we are dealing with, the final result for the 
gravitational wave spectra would be the same, since there are no lack of information 
by matching only a part of the metric perturbation, the reduced amplitude. The 
remaining part would give us redundant information. 



But in the case of theories with dilaton fields, Brans-Dicke frames and in par- 
ticular, string cosmology, the remaining part in the total amplitude perturbation 
contains additional information with respect to the reduced one. It contains an 
added dependence on the dilaton field, as can be seen from eq. ([4.16|) . The be- 
haviour of this part through the transition is not the same as the scale factor, since 
there is no continuity in the first derivative of the dilaton field. Thus, the power 
spectra computed by matching only the reduced amplitude perturbations is different 
from that obtained by matching the full metric perturbation. 



In summary, the physical quantity to be studied is the full metric perturbation 
amplitude. Meanwhile, the reduced amplitude perturbation can be considered as a 
sometimes convenient mathematical reduction of the problem. 



In the inflationary stage, we have the same equations for the reduced amplitude 
perturbation ( [4. 131 ) and ( [4. 141 ) that in the precedent section. Therefore, we will have 
the same expression for the reduced amplitude perturbation: 

^/ = U{m - ftYHu^^^ ik{r]i - f])) (5.54) 

with 

-4^^ 4) (-^) 



The reduced -0 and the total tensorial metric perturbation are related by eq. (|4.15| ) . 
Instead of ^, we will match the metric perturbation amplitude $: 

d-l 

$ = 2 e2 (5.56) 
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such that the metric perturbation is written as: 

V = e^,(k)$(?7,k)e^^-^ (5.57) 
From eqs.( p.6|) , (|3l5| ) and ( |5.47|) , we have: 

$^ = Uai^e^ {r]j - 7^)"^(^)+57^^W (k{r]j - fj)) (5.58) 



In contrast to the precedent cases, we will take in the following the negative sign in 
the parameter i> ( 5.55| ): 

because with this definition, we can write the amplitude perturbation as: 

(tj = Afai^e^ {r]j - nfUy^^^ {k{7]i - fj)) (5.60) 

where A/" is a normalization constant. Notice the parameter u in the conformal time 
dependence. 

On the other hand, in the radiation dominated stage, the reduced amplitude 
perturbation is obtained in a totally equivalent way to the precedent section. We 
recover for this reduced amplitude eq.( ^.53|) : 

since we have again the parameter = ±|. 

The tensorial metric perturbation h^u is given by eq. (|5.57|) with the amplitude 
$ given by: 

= tM^e^aii-^f]-<'-^) (ae-''^^*?-''^) + Pe""^"^-^'^) (5.62) 
where r = In shorter form, introducing the notation: 

, d+l </>/ 

j^ = J\fai — e- (5.63) 
g = iJ- M air^e^ , (5.64) 

V TT 



we have: 



= J^(nr - fiyn,^'^ (kirjj - m (5.65) 
= gk-'^fj'^ [ae'"'^''"'"'^ + /3e*^(^-^'i)) (5.66) 
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where v is given by eq.( |5.55| ). 



As indicated above, in this case we will compute the Bogoliubov's Coefficients 
by matching the total amplitude of tensorial metric perturbation $. Thus, we have 
to solve the equations: 



which yield: 



dr] 



Vi) 



V = Vl 



^iiiv = Vi) 
d^ii 



dr] 



n = ni 



gk'^rii'^{p{l - ikrii) + a(l + ikrn)) 



(5.67) 
(5.68) 

(5.69) 
(5.70) 



By solving this system, we have: 



^-f + t]n.^'\x) + n'J'\x) 



X X 



nJ'\x) + n'J'\x) 



(5.71) 
(5.72) 



where X is given by eq. ( |5.19| ) 



Finally, after computation and using the properties of Hankel functions and 
their wronskian (See and pSf), we have the following expressions for the square 
moduli: 



a 



2u+2 



q\'^'' k 



,g\^ " Vr/ 4 

+ n.^,^'\x)n._,^'\x) 



I 



X2 



nj^^\x)nj^'\x) 



(2), 



(5.73) 



/5 



^1 



4 



+ 7^,_i«(X)7^,_i(2)(X) 



2 ZL 2 zL 



(5.74) 
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From expressions ( |5.73| ) and ( [5 .741 ) is easy to see: 



a 



13 



Q 



2u+2 



q\^'' k 



n k r]i 



(5.75) 



which by applying the normahzation condition of Bogohubov's Coefficients eg. .311) , 
gives: 

^r^,-+^f^)""^ = 1 (5.76) 



TV 



Now, we can translate this condition over the usual normalization constants J\f 
and Ad. From the definitions (|5.63|) and ( p.64 ): 

_ lATf TV e^' 

From the matching properties of the scale factor in conformal time eq. 

aij 

And from eqs. (|313|) and (CT) : 



r+q 



e'Pll 



^ -2d -2dr 
ail Vl 



With eqs. (|5.77|) and (|5.78|) and the u value eq.( |5.59|) we obtain: 



JV \ I q\ l-2u+r{l-d,) 



Since r = eqs.( p.76D and ( |5.79| ) yield: 



M 



Vl 



2-r{l-d) 



we have: 
(5.77) 

(5.78) 

(5.79) 



(5.80) 



That is, the normalization factors A/" and Ai differ only in a phase. We can 
eliminate the normalization constants in the expressions for the Bogohubov's Coef- 
ficients. From eq.( |5.76 ), the coefficient multiplying the square brackets expression 
in eqs.( ^.73"D and ( ^.74| ) is exactly ^^kr]i. Then, the final expressions for the a and 
l3 coefficients are: 



|2 ^ V 

a \ = — X 



1 + ^ 

V 

n,^,^'\x)n,^,('\x) 



+ ^ \ n,^'\x)n,^'\x) + 



- 2j^nJ'\x)n.^^^'\x)~^^^ + -^ 



(5.81) 
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(3 



X 



X2 



+ n,^i^'\x)n.^^'^'\x) 



(2)/ 



X 



ttX^ ttX 



(5.82) 



5.3.1 The Power Spectrum in the Full Dilaton Case 



Using eq. (|4.29| ) and ( ^.82| ), we have for the power spectrum the next expression: 



P{uo)du 



^TTC^ S 
\2 



\ 



+ ^7^.-i«(c.5)7^.-i(^)(a;5) + 
- 2^7i.«(u;5)7^._/^)(c.5)-^A ^ 



7r2 7r^-^2 



(?)■ 



(5.83) 



The fraction of critical energy density takes the expression: 

, 2 



hG 



1 + ^^|H.'"('^S)H„'«(u;S) + 



+ 



2^^7^.«(a;5)7^._/^)(a;5)-^A ^ 



(5.84) 



We compute for the three-dimensional case the values of the coefficients in front 
of the square brackets in eqs. (|5.83|) and ( ^.84| ). These values differ from our previous 

computations on ( ^.39| ) and ( ^.46| ) in a factor (^^;zx^ • Thus, in the full dilaton case, 
we have: 



87rc3 S 
h G 



4.008 10 



-55 erg seg" 
cm^ 



6.220 10"^^ seg-^ 



(5.85) 
(5.86) 
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In the Tables (^) and (^), we show the different power spectra and fraction of 
critical energy density for a single cosmological scale factor, but in the three different 
roles for the dilaton. As can be seen, there is a notorious difference between the 
three different treatments of the dilaton role for a single cosmological scale factor 
description. First, the parameter u, characterizing the inflationary expansion and 
the power spectrum, is different in each case. The sign of this parameter is only 
ensured in the full dilaton treatment, by relating it with the form of the metric am- 
plitude perturbation. Second one, the coefficients in front of the square brackets are 
different. Within the reduced amplitude treatment (no dilaton and partial dilaton 
cases), these coefficients appear related with the quantity ('^^ l), meanwhile for 
the total metric amplitude hi^^ matching (full dilaton case), the coefficients have the 
expression ^2^. 

It must be observed that the No Dilaton Case is coherently described also by 
the expressions corresponding to the Full Dilaton Case. In fact, in the first one, the 
next relation holds: 

^ 2 r 

(this relation is not satisfied by the partial dilaton role case, which must be described 
by the specific expressions above indicated). 

These observations enable us to work with general expressions for the stochastic 
background of gravitational waves, given by the corresponding ones of the full dilaton 
case eqs. (|5.83|) and ( |5.84|) . The partial dilaton and no dilaton cases can be simply 



obtained from these expressions by substituting the factors ^ by [i' — |j with the 
proper value of parameter u for each one. 
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Case 



par. inflationary v 



gral. 



ci = 3 



{y-\) 



Power Spectrum 
P{uj)duj 



No 
Dilaton 



87rc3 



du 



Pctrtial 
Dilaton 



2 u)S 



Full 
Dilaton 



(i - </¥) 



n , ,2 ill 



■(io; 



SY 
{if 



4 u}S 



Table 3: The Dilaton Role in the Power Spectrum in String Cosmology 

Comparative Table of the Gravitational Waves Power Spectrum obtained in our String 
Cosmology Model. Differences have oughted to the weight of the dilaton in each case. 
The differences between the three cases turn out in the value of the u parameter and the 
coefficients in the power spectrum expresssions. 
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Case 



par. inflationary v 



gral. 



d = 3 



iy-V) 



Contribution to Energy Density 



No 
Dilaton 



h G 



Pcirtial 
Dilaton 



± 



1 + 



2, 

4:UJS 
I,, \\ 



Full 
Dilaton 



n G 



3Ho'^c^ S 



CO 



1 + 



sy 

w 



Table 4: Contribution to Energy Density with different Dilaton Role in String 
Cosmology 

Comparative Table of the Fraction of the Energy Density ^gw contained in Gravitational 
Wave Backgrounds produced in our String Cosmology Model. Again, the differences at 
the value of v parameter and have oughted to the weight of the dilaton role in each case. 
Signatures of this are the value of the v parameter and the coefficients in the power 
spectrum expresssions. 
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6 Asymptotic Behaviours of the Gravitational Wave 
Background 



We have computed the analytical expressions for P{uj)duj and Qqw at high 
and low frequencies. Both regimes correspond to very large {X >> 1) and very 
small {X << 1) arguments of Hankel functions, respectively. We consider the 
three regimes: low frequencies for uj << {X << 1), medium frequencies for 
a; ~ 1 (exact analytical expressions ( p.38[ ),( p^^ ),( |5.8^^ ),( p.84D ) and high frequencies 



for a; >> {X » 1) where is the characteristic frequency (eq.( ^.36|) ) with 
order of Mhz for the scale factor evolution considered. 



The exact analytical form of the power spectrum, eqs.( ^.38D and ( ^.83|) requires 



a careful algebraic treatment in order to compute their asymptotic behaviours. At- 
tention must be put both in the range of validity as on the order of expansions used. 
Computation is quite long, thus we omit details here and give only the leader orders. 



6.1 High frequencies 

The High Frequency Asymptotic expressions for the No Dilaton (ND) and 
Partial Dilaton {PD) Cases are formally equivalent. The leader orders for the power 
spectrum and contribution to energy density can be written respectively as: 

P{u)duU^^^ - '-^u-' + 0{u-') (6.1) 

^o.U,D - + Oi.-^) (6.2) 

The Full Dilaton {FD) Case presents different expressions for the High Asymp- 
totic Behaviour: 



2hG (^) 



2 



^Gwl^n - TT^V^^^' ( ^ - 1 ) + Oiuj') (6.4) 



Notice that leader order terms have fixed dependences, independent of number 
of spatial dimensions or parameters of the solutions. In fact, P{uj)du in both No 
Dilaton and Partial Dilaton cases present a convergent behaviour at high frequencies 
as ~ u~^. In contrast, the Full Dilaton case is divergent in this regime as ~ a;^. 
Vtaw in the No Dilaton and Partial Dilaton cases reachs asymptotically a constant 
value. In the Full Dilaton case VLqw diverges as u"^. 
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6.2 Low frequencies 



Opposite to the high frequencies case, leader orders for low frequencies depend 
explicitly from the parameter v. As a consequence, each case presents a particular 
expression for their asymptotic behaviours. To leader order uniquely, the No Dilaton 
case and the Partial Dilaton case have the same formal expressions: 



2 



P{u)du\^^^p^ ~ 3 ^ du—TiuYiuSr''' + OiiuSY) (6.5) 



^GwU,PD - ^Tiu)\uSf~"' + 0{{uSr) (6.6) 



hG [^-i) 2^^ 



In the three-dimensional case, = | for the (ND) case and z/ = | for the 
(PD) case. We have at very low frequencies P{u)duj\j^^ ~ and ^gw\nd ~ 
P{u)du\pjj ~ U3 and flcwlpD ~ 

In the Full Dilaton Case the low frequencies behaviour is: 

Pi-)''-\.o ~ ^|1-<^^W^ (1 - y + OH.Sr) (6,7) 
-l^M ^^(l - -Xi-Sr^-^OH^Sf, (6.8) 



The leader order has the same behaviour P(uj)du ~ ojs and ilcw ~ tus as in 
the Partial Dilaton Case, but the coefficient is different. 

The leader order at very asymptotic behaviours for all cases is summarized in 
table d). 



Case 


P{uj)duj 








Low Prq. 


High Frq. 


Low Frq. 


High Prq. 


String 


No Dilaton 


i 




4 


cte. 


Driven 


partial Dil. 


~ 




8 

~ U3 


cte. 




Full Dil. 


b 


~ u 


a 

~ a;3 





Table 5: Asymptotic Behaviours of Gravitational Wave Power Spectrum and Con- 
tribution to Energy Density. 

The dominant dependence for very low and very high frequences are summarized. The 
values are refered to the three spatial dimensional case. 
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7 Discussion and Conclusions 



We have studied the production of a primordial stochastic gravitational wave 
background in a cosmological model fully extracted in the context of selfconsistent 
string cosmology. As seen in Section II, suitable descriptions for inflationary, ra- 
diation dominated and matter dominated stages can be extracted from effective 
string theory in which strings propagating in the curved backgrounds are the clas- 
sical sources of matter and drive selfconsistently the cosmological background. The 
equation of state provided by the dynamics of string matter naturally generates in 
its evolution the background in the three stages. 

No exact description of the transitions among stages can be extracted in the 
framework of this effective treatment. Higher order corrections in low energy effec- 
tive equations or more detailed study on the evolution of the gas of strings would 
be neccessaries in order to handle such subject. 



This lack of full understanding of the transition dynamics is an usual limita- 
tion when primordial metric perturbations are computed. The main consequence 
is the loss of predictibility. Usual gravitational wave results depend on parameters 
of the cosmological model which are not always well known whose 



physical meaning must be recovered after computation (||16|), or where the not con- 
sidered transition dynamics must be recalled in order to overcome divergences (0], 
Usual arguments dealing with horizon exit and reentry require a 

suitable cosmological description in order to be translated univocally into current 
measurement abilities. 



We faced this problem before to compute the metric perturbations. We elabo- 
rate a cosmological minimal model of evolution for the scale factor with satisfactory 
sudden and continuous transitions, provided it makes use of particular descriptive 
cosmic time variables at each stage (see eqs.( p.2|) ). In this way, we merge our lack 
of knowledge about real transitions by modelizing them in a descriptive scale factor 
expressed in suitable cosmic time variables. This description is linked with the min- 
imal observational Universe information (standard values for transition times and 
scale factor ratii reached in each stage). It provides an equivalent evolution for the 
scale factor, but evades the undetermined region of transitions without introduction 
of free parameters or loss of predictibility. Unless the details imprinted by dynamics 
transitions, in our case modelized as sudden and continuous, our computation is 
expected equivalent to one made on a full physical model running on cosmic time. 

We have also imposed continuity and smothness on the conformal time descrip- 
tion, by constructing this on the descriptive variables (see eqs.( |3.5| )). Continuity 
of scale factor both in cosmic time and conformal time description is a feature not 
always present in literature (see |^). In our case, the variables constructed en- 
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able us to maintain the linking with the observational Universe information. This 
model seems us more appropiated in order to perform primordial metric perturba- 
tions computations, since no effects oughted to discontinuity in the scale factor of 
the metric are introduced. Not so perfect is the treatment of the dilaton field, whose 
evolution can be treated solely continuous but not smooth. Further study would be 
needed in order to valutate this effect. 



The only gravitational waves we have computed are those produced as ampli- 
fication by the evolution of background of tensorial metric perturbations. We do 
not account here another possible mechanisms of generation neither amplification 
of another metric perturbations. 



The variable X in the power spectrum and the proper frequency u are related 
in a way totally determined by the description of the cosmological scale factor evo- 
lution, see eg. ( ^.34 ). The factor relating them depends of the expansion ratii, the 
exit time of inflationary epoch and the coefficients of expansion inflationary and 
current epochs. Being all them fixed in our cosmological background linked to the 
observational Universe, no free parameters are introduced at this level. No one of 
remaining unknown parameters, like global scale factor A//, appears on the results 
of our computation. Differently from almost all string cosmology computations in 
literature, firm predictions on precise frequencies ranges can be extracted in our 
case. 



In this way, we have computed exact, fully predictive and free-parameter ex- 
pressions for the power spectrum P{uj)du and contribution to energy density VLqw 
of the primordial gravitational waves background generated in the transition among 
the inflationary and radiation dominated stages. We have not considered the contri- 
bution given by the radiation dominated-matter dominated transition in this study. 
The gravitational wave contribution due to this transition is expected to be ne- 
glectelly small, as compared to the first transition. It would be expected the second 
transition having a role only on the low frequencies regime [Q, not so important in 
anycase for our results. 



7.1 The Dilaton Role 



We have obtained drastic differences in the stochastic background of gravita- 
tional waves produced in the same scale factor evolution by considering differents 
degrees of the role played by the dilaton. The simplest case, without account of 
the effect of the dilaton neither on the perturbation equation nor on the amplitude 
perturbation. The second partial account, with the proper perturbation 

equation but still matching the reduced amplitude perturbation. The lastest case, 
a full account by working with the total tensorial amplitude perturbation and per- 
turbation equation. 
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The background of gravitational waves comes characterized in their shape by the 
parameter z/, which depends of the inflationary description, the inflation-radiation 
dominated transition (the only on which we are computing gravitational wave pro- 
duction) and the role played by the dilaton. The expressions for v have been found 
in the three cases (eqs.( ^.4|) , (|5.49|) and (|5.59|) ). We obtain an exact expression for 
the power spectrum and energy density contribution (eqs. ( pT55D and ( |5.45D ) in terms 
of Hankel functions of order z/, formally equal in the No Dilaton and partial Dilaton 
cases. Differences among them are oughted to v differences. The formal expressions 
in the full dilaton case are different (eqs. (|5.83|) ,( ^.84| )) both in parameter v as in 
coefficients involved. 



The low frequency and high frequency asymptotic regimes have been discussed. 
In the No Dilaton case, asymptotic behaviours for power spectrum are both vanish- 
ing at low and high frequencies as u;3 and uo~^ respectively. This gives a gravitational 
wave contribution to the energy density asymptotically constant at high frequencies 
of magnitude O^gw ~ 10~^^. There is a slope change that produces a peak in the 
power spectra around a characteristic frequency totally determined by the model of 
uo^ ~ 1.48 Mhz. 

The Partial Dilaton case introduces the effect of the dilaton exclusively in the 
tensorial perturbation equation (which is not longer equivalent to the massless real 
scalar field propagation equation but not on the perturbation itself. The 

general characteristic are very similar to the previous case. Both asymptotic regimes 
for P{uj)duj vanish again, but with dependences us and u~^. The peak appears 
around the same characteristic frequency, with value one order of magnitude lower 
than in the No Dilaton well as the asymptotic constant contribution to 

energy density. 

In contrast, when the full dilaton role is accounted, general characteristics as 
well as orders of magnitude of the spectrum are drastically modified. It has similar 
values for the frequencies below the Mhz, with power spectrum vanishing again as 
LOS. For high frequencies, in opposition to the former cases, both P{u)du and flaw 
are increasing at high frequencies. For P{uj)du, an asymptotic divergent behaviour 
proportional to u is found. It gives values much higher than the no dilaton and par- 
tial dilaton cases. The contribution to ^gw is equally divergent at high frequencies 
as w^. The change of slope is less visible and no clear peaks are found. The transition 
from the low frequency to the high frequency regime is slower than in the previous 
case and the full analytical expressions are needed on a wider range 10^ ~ lO^Hz. 
Comparative tables of these summarized results can be found in tables (|^), (|) and 
(1). See also figs.(P and (g). 

Existence of an upper cutoff must be studied, whereas the Full Dilaton Case 
studied could be in contradiction with observational bounds as the current value of 
total energy density in critical units ~ 1. In that case, an end-point not predicted 
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by the current minimal model considerations could be introduced in the spectrum as 
made in the literature [|18|. Divergent high frequencies behaviour and introduction 
of an upper cutoff is an usual feature in the string cosmology contexts. 



7.2 String and No-String Cosmologies 



Among the spectra computed in string cosmology contexts, it must be distin- 
guished between those computed in Brans-Dicke frames (that we compare with our 
Full Dilaton Case) and those computed following an usual quantum field theory way, 
that is, as our No Dilaton Case. 

The shape of the spectra computed in string cosmology contexts are very sim- 
ilar. The principal features, as slope changings, are signal of the number of stages 
or transitions considered in the scale factor evolution. All the known cases, co- 
herently treated in Brans-Dicke frames, presents an increasing dependence at high 
frequencies. 



Notice that we are not handling with a "Pre-Big Bang" scenario |14|. Our String 



Driven Cosmological Background runs on positive values for the proper cosmic time 
t. It have been said |T^ that Pre Big-Bang scenario predicts a power spectrum with 
a peak and an end-point. Our analysis enable us to identify these features with the 
Brans-Dicke framework where the low effective energy treatment is made, without 
necessity of a "Pre Big Bang" phase. In fact we have seen as the same Universe 
evolution and transitions have given different power spectra. In the full dilaton 
treatment, the transition among the inflationary inverse power stage and the radia- 
tion dominated stage gives an always increasing function for the number of particles 
causes the high frequency range to be asymptotically divergent, feature 
that disappears completely in the No Dilaton or in the Partial dilaton treatments. 



Similarly, a cosmological model involving successively a named Dilaton Driven 
stage plus an "String Phase" at nearly constant curvature, and a radiation dom- 
inated stage produces also a spectrum increasing with frequency ([0). As made 
in the literature (|]l8l), an upper cutoff is placed on this increasing spectra by con- 
sidering the frequency making | (3{u>) |^ = 1 as the maximum one produced. This 
frequency is the one where one graviton per space phase unit volume and polar- 
ization is produced, and after which exponential suppression is supposed. These 
arguments leds to identification of the so called "end-point" that acts as "peak" 
of the spectrum oughted to the increasing power with frequency. The exact values 
at current time are functions of the undetermined parameters in the cosmological 
model there treated 
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We conclude that whatever gravitational wave computation on inflationary 
stages of the type extracted in string cosmology, coherently made in the Brans- 
Dicke frame, must give an increasing spectrum. The peaks are produced by slope 
change and they are signal of the transitions in the dynamics of background evo- 
lution. We consider the Pre-Big Bang scenario do not predicts a peak, but it is 
supposed by defining a ui such that |/3(a;i)p = 1. This proper frequency is com- 
puted at the beginning of radiation dominated stage, when the wave reenters the 
horizon and it must suffer a redshift at current time, expressable as function of 
unknown parameters of a "string phase". It acts as end-point because waves with 
uj > uji are supposed exponentially suppressed. Since the spectrum was increasing 
with frequency, the same frequency constitutes a maximun (peak). 

If we use the same argument in order to fix an upper limit, our spectrum must 
be cutted at frequency uJmax ~ 3.85 MHz where the power spectrum will have a 
value aroimd P{uj)dLO ~ 5.68 lO^^^fff and now ~ 3.40 IQ-^^pc- These are the 
same order of magnitudes of the peaks atteinted on the No Dilaton and partial 
Dilaton Cases. No confiict with observational constraints look possible for such 
weak signals predicted. But this argument could be too naive, since in the practical 
way is equivalent to cancelate from the spectra the features introduced by the full 
dilaton role. 

For instance, if we want to relaxe this condition and we ask us what would 
be the gravitational wave with maximum possible frequency, we can give at least a 
coherency condition. Now, let us ask if it is possible to produce a wave with period 
lower than the Planck time, that means with wavelenght lower than the Planck 
lenght scale. If that feature is possible, it is clear that our simplified quantum the- 
ory formalism would be unappropiated in order to compute such produced particles. 
We can fix the upper cutoff in the corresponding cup ~ W^Hz by considering that 
beyond it the power spectrum of gravitational waves must be, at least, different 
from elsewhere hitohere computed. At those values, both the No Dilaton and par- 
tial Dilaton case predicts completely negligible production, but not is the case for 
Pull Dilaton Case. Its asymptotic divergent behaviour predicts enormous values for 
P{uj)du} and ^gw- 

Obviously, we do not mean that these are the right predictions of string cosmol- 
ogy. Without doubt, the model here treated is excessively simple in order to handle 
properly the end-point of spectrum. But usual assumptions on upper exponential 
suppression could be more careful revisited, since a too naive treatment in this point 
can mask important features of predictions. 

In relation with no-string infiationary cosmologies, the so called standard infia- 
tion is usually intended as a De Sitter stage. Notice that string driven inflationary 
stage describes an evolution with inverse power dependence. It must not be confused 
with usual power law, although our model can be said superinflationary too. 
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There is a radical difference among the string cosmology spectra and those 
obtained with an exponential inffationary expansion [Q. The divergence at low 
frequencies that it supposes is not found in our String Driven Cosmological Back- 
ground The high frequencies behaviour are compatible with the No Dilaton Case 



(|5.38|) , computed in a totally equivalent way. 



If comparaison is made with the Full Dilaton Case, a totally different behaviour 
is found with respect to the obtained in De Sitter case. Notice that comparaison 
among string comology inflationary models and standard inflation means confrontate 
power and inverse power type laws with De Sitter exponential inflation, since until 
the moment no De Sitter type expansion have been coherently obtained in String 
Cosmology. This difference in the inflationary scale factor dynamics, together with 
the appropiated teatment of metric perturbations in each framework are the main 
causes of differences among the gravitational wave power spectra in both cases. 

No explicit dependence on the beginning of the inflationary stage has been 
found on the gravitational wave computation. In anycase, further study must 
be done in order to determine the influence on the power spectra if earliers in- 
flation stages are considered. Gravitational wave production in frameworks with 
multistages inflation can be found in literature, both in the most classical cases 
(p6[,P|,0) meanwhile in the string cosmology contexts there are various multistage 
models 



In summary, there are a variety of effects affecting the gravitational wave power 
spectra and justifiying the differences among our results and another string and no 
string cosmologies. The most notorious difference is oughted to the existence of a 
dilaton field and a metric evolving in a Brans-Dicke frame. Differences are found 
on the equation for the tensorial perturbation, and the amplitude of perturbation 
depends itself on the dilaton field. Our spectra here computed have showed clearly 
the successive effect obtained on the gravitational wave production by introducing 
the dilaton effect on the perturbation equation (partial Dilaton case) and on the 
perturbation itself (Full Dilaton case). From an unique evolution of the scale factor, 
we have seen as the dilaton role on the spectra consists in increasing the high fre- 
quencies regime, leading to a divergent behaviour. This is in agreement with almost 
string cosmology, Brans-Dicke computations. In our cases without full account of 
the dilaton role, we have spectra decreasing at higher frequencies, as such obtained 
in the literature |Q in the more standard frameworks. 

There are many points that desserve further study in order to be clarified, either 
from the point of view of string cosmological backgrounds and transitions dynam- 
ics, either from those of gravitational wave computations in the string appropiated 
frameworks. Better treatments than here applied could take place in every phase of 
the problem, advising us against consider this procedure or its results as definitives. 
In anycase, we have proved here some aspects of the way predictions on observational 
consequences can be extracted from String Cosmology. 
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